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I. INTRODUCTION 



In this paper the evolution of structural relaxation as observed upon cooling the van- 
der-Waals liquid propylene carbonate (PC) from above the melting temperature (T^ = 
218 K) to the glass-transition temperature {Tg = 160 K) will be analyzed. It will be shown 
that the spectra, as measured within the four-decade frequency window below 800 GHz by 
depolarized-light-scattering, by dielectric-loss, and by neutron-scattering spectroscopy can 
be quantitatively described by the solutions of a two-component schematic model of the 
mode-coupling theory (MCT), where the drift of the various spectral features over several 
orders of magnitude due to temperature changes can be fitted by smooth variations of 
the model parameters. The results of the data fits will be used to demonstrate in detail 
which features can be explained by the universal /3-relaxation-scaling laws of the asymptotic 
MCT-bifurcation dynamics, and which are caused by either preasymptotic corrections to 
this scaling or by crossover phenomena to microscopic oscillatory motion. 

Glassy PC spectra within the full GHz window have first been studied by Du et al. 
Ill] using depolarized-light-scattering spectroscopy. It was shown that the data can be in- 
terpreted with the universal laws predicted by MCT. In its basic version, which is also 
referred to as the idealized MCT, this theory implies an ideal liquid-glass transition at 
a characteristic temperature T^. In an extended version, Tc marks a crossover from the 
high-temperature regime, where the dynamics is dominated by non-linear-interaction effects 
between density fluctuations, to a low-temperature regime, where the dynamics deals with 
activated-hopping transport in an effectively frozen system. For temperatures T near Tc, 
the MCT equations can be solved by asymptotic expansions for the so-called /3-relaxation 
regime. This results in formulas for universal features of the MCT dynamics as reflected in 
the appearance of dynamical scaling laws, power-law-decay processes, and in algebraically 
diverging time scales. The different anomalous exponents and also the /3-relaxation-master 
functions are determined by a system-dependent number which is called the exponent pa- 
rameter A [^. The data analysis of Ref. |]1| suggested T^. ^ 187 K and A ~ 0.78. Relaxation 
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curves measured for PC within the pico-second window in solvation-dynamics studies 
and dielectric-loss spectra determined within the GHz window have also been ana- 

lyzed with the MCT-scaling-law formulas using parameters Tc and A consistent within the 
experimental uncertainties with the values cited above. The critical temperature for PC has 
first been determined to ~ 180 K |^ by interpreting the a-relaxation time for density 
fluctuations measured by neutron-scattering spectroscopy with the MCT-power-law predic- 
tion for this quantity. A similar analysis of the viscosity |l|,0] suggests a value of Tc near 
190 K. The effective Debye- Waller factor for the elastic modulus has been measured for PC 
by Brillouin-scattering spectroscopy ||^ . Interpreting this quantity with the asymptotic for- 
mula of the idealized MCT, a critical temperature considerably higher than 190 K has been 
suggested. However, since the data interpretation is not compelling |l[], this finding cannot 
be considered to be a falsification of the Tc ~ 187 K result. Thus one could conclude, that 
MCT describes some essential features of the glassy dynamics of PC qualitatively correct, a 
statement which also holds for a series of other glass- forming systems |0. 

In order to arrive at a more stringent assessment of MCT, Wuttke et al. have re- 



examined the above cited PC data for T > Tc. In addition, they have studied incoherent- 
neutron-scattering spectra S{q, u) for a two-decade window in frequency u and for wave 
vectors q between 0.7 and 2.3 A ^. The data exhibited the predicted factorization in a 
g-dependent but cj-independent amplitude hg, and a g-independent term describing the 
frequency and temperature variation: S{q,u!) oc hqx" iyj) / u . The susceptibility spectrum 
x"iy>) showed the subtle dependence on uj and on (T — Tc) predicted by the MCT-scaling 
laws for the /^-process, provided Tc ~ 182 K and A ~ 0.72 was chosen. These parameters 
are marginally compatible with the values found in the above cited earlier work on PC. The 
depolarized-light-scattering spectra have been remeasured within the /3-relaxation window 
for T > Tc- The spectrometer used in Ref. |T0| incorporated several improvements over 



the one used in the original study resulting in improved signal-to-noise ratios. Further- 
more, the use of a narrow-band interference filter eliminated the possibility of higher-order 
transmission effects which have recently been recognized as a potential source of artifacts 



rT| , [T2[] - But the new spectra agree with the old ones within the error bars of the latter. 
The remeasured spectra could be fitted convincingly with the universal asymptotic results 
using the newly found values for Tc and A. It was shown in addition that also the solvation- 
dynamics results and the dielectric-loss spectra could be fitted within the same 

frame using the new values for and A. Actually, the new fit to the dielectric-loss data |T^ 



is more convincing than the original one 0-^], since the fit interval expands with decreasing 
(T — Tc), as requested by MCT. The size of the (T — Tc) interval and the window for the 
frequency where leading-order-asymptotic results describe the MCT-bifurcation dynamics. 



depend on the probing variable |jT3|JTj]. It was assumed in [jTO[ that the range of validity of 
the asymptotic analysis is smaller for the dielectric-loss spectra than for the light-scattering 
spectra. It also had to be anticipated that preasymptotic corrections can account for a 35% 
offset of the /3-relaxation-time scale of the neutron-scattering data relative to the one for the 
light-scattering data. 

To corroborate the cited MCT interpretations of glassy PC spectra, the previous work 
shall in this paper be extended in three directions. First, the a-relaxation peaks will be 
included in the analysis, so that the low-frequency limit for the fit interval can be decreased 
to 1 GHz or lower. Thereby the crossover from a- to /^-relaxation and the non-universal a- 
peak shapes can be described as well. Second, the crossover from relaxation to vibrational 
dynamics will be included in the analysis, so that the high-frequency limit for the fit interval 
can be increased by about a factor of four. Third, an extended form of the MCT instead 
of the idealized one will be used, so that the spectra for depolarized-light scattering and 
dielectric loss for > 1 GHz can be described also for temperatures below Tc. The specified 
goals will be achieved by studying the full solutions of an MCT model. 

The paper is organized as follows: In Sec. y, the basic formulas for the schematic model 
to be used will be summarized, and then (Sec. fT|) the experimental data sets are fitted using 
this model with smoothly drifting parameters. After a short introduction to the necessary 
equations for the asymptotic analysis for the model (Sec. fV|), the /^-scaling laws are tested 
against the data in Sec. 0. In Sec. 0, it will be shown that for the studied model a properly 



defined dielectric modulus is more suited for a description by scaling laws than the dielectric 
function. Section |V11| presents some conclusions. 



II. A SCHEMATIC MODE-COUPLING-THEORY MODEL 

The idealized MCT is based on closed equations of motion for the auto-correlation func- 
tions of the density fluctuations (pqit), which are positive definite functions of time t, de- 
pending on the wave vector modulus q |T5|. The extended MCT also includes couplings of 
the density correlators to the auto-correlation functions for the currents |TB[. The gen- 



eral equation of motion expresses the density correlator in terms of relaxation kernels. It is 
formulated most transparently with Laplace-transformed quantities. For the latter, the con- 
vention F{z) = i Jl^ exp{izt)F{t) dt with complex frequency z, and F{uj) = F'{uj) + iF"{uj) 
ioT z = oj + will be used. 

^'M-^'W- .+Mr'("ln;^.(.r '''' 

Here, Vtq denotes a characteristic frequency given by the thermal velocity v and the static 
structure factor Sq\ fi^ = q^v"^ / Sq. The general current-flow kernel Cq{z) describes density- 
fluctuation decay via two parallel channels. Phonon-assisted hopping is given by Nq{z). 
The relaxation due to nonlinear interactions of density fluctuations is described by a force- 
fluctuation kernel which consists of a sum of a regular term M^'^^{z) and a mode-coupling 
term mq{z). The former deals with normal-liquid dynamics, and the latter with the slow 
motion caused by the cage effect. It is obtained as a polynomial J^q of the density correlators 
<Pq{t): 

mq{t) = Tq \4>q{t)\ . (ic) 

The coefficients of the polynomial are non-negative; they are given by the equilibrium struc- 
ture and hence depend smoothly on external control parameters like temperature T. Sys- 
tematic studies of the kernels Nq{z) and M^'^^{z) are not available. The theory shall be 
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simplified by Markov approximations of tliese quantities: M^'^^{z) = iuq, Nq{z) = iAq. 
The friction constants z/g > and flopping coefficients > sliaff be treated as model 
parameters, which depend smoothly on T. 

Equations (0) can exhibit bifurcation singularities. Generically, if as a single control 
parameter the temperature is considered, the singularity occurs for a critical temperature 
Tc if all hopping coefficients Aq vanish. If some 7^ 0, the singularity is avoided. However, 
for small Aq and small \T — Tc\ the singularity causes an anomalous dynamics: the glassy 
dynamics studied by MCT. At the singularity the correlators do not decay to zero but to a 
positive value fq, which is called the plateau. It is approached by an algebraic decay law, 
called the critical decay, which is specified by an anomalous exponent a, < a < 1/2: 

Ut) - r, = hqit/to)-' + O (t-2») ; 

T = T„ A, = . (2) 

The quantity /ig > is called the critical amplitude, and it can be determined from the 
mode-coupling functional J^q for T = T^. The time scale to is determined by the transient 
dynamics for T = T^. For Aq = and small but negative (Tc — T), the correlator falls below 
the plateau fq according to the von Schweidler law (f)q{t) — oc —t^ + O (t'^''^ , characterized 
by a second anomalous exponent b, < b < 1. From J-'q for T = T^, one can calculate the 
above mentioned exponent parameter A, < A < 1/2, which determines the critical exponent 
a and the von Schweidler exponent b via r(l — a)^/r(l — 2a) = A = r(l + 6)^/r(l + 26). 
In the so-called /3-relaxation window, implicitly defined by \(pq{t) — fq \ <^ 1, MCT predicts 
that the dynamics is in leading order controlled by merely two smooth functions of T: the 
separation parameter a and the hopping parameter 6. The former is determined by J-'q, and 
its zero defines the crossover temperature T^. a = C{Tc — T)/Tc, + O ((T — Tc)^). The latter 
obeys 5 > 0; generically, 5 vanishes only if Ag = for all q. The shape of the correlation 
functions in the asymptotic regime of the /3-relaxation window is fully determined by the 
exponent parameter A; as can be inferred from Ref. |0] and the original papers cited therein. 
Testing the relevance of MCT by comparing the leading-order results for the /3-relaxation 
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with data is however hampered by a great difficulty. Without detailled microscopic calcu- 
lations one cannot determine the size of the corrections to the asymptotic formulas, and 
therefore their range of validity is not known. In addition, the optimal choice of A, fixing 
the shape of the log x"-versus-log oj graph is tedious to decide upon and might well depend 
on the choice of the fit interval. The difficulty of fixing A from a /3-relaxation study alone 



was demonstrated recently for the hard-sphere system fl^. A set of density correlators 
calculated for various wave vectors and packing fractions was considered. A ffi to them with 
the asymptotic predictions for a significantly wrong A was by a standard fitting procedure 
not distinguishable from the correct ffis within typical experimental windows. 

A different route for data interpretation is based on comparison of the measured spectra 
with the complete solutions obtained from schematic MCT models. This procedure was 



studied ffist by Alba-Simionesco et al. [18-20|. Schematic models are truncations of the 



complete set of Eqs. ([^) to a set dealing with a small number of correlators only. Thus 
the mathematical complexity of the problem is reduced considerably. Alas, the connection 
of the mode-coupling-functional coefficients with the microscopic structure gets lost; the 
coefficients are to be treated as ffi parameters. The main advantage of this approach is that 
one does not rely on the applicability of asymptotic formulas; one is sure that all results on 
crossover phenomena and preasymptotic corrections are logically consistent with the MCT. 

The simplest schematic model deals with a single correlator only, which shall be denoted 
by (pit). The ffist MCT equation is equivalent to Eqs. (0,b) with q indices dropped: 

'^^^^ " z + i/\-n^/[z + iv + n^m{z)\ ■ ^^^^ 

For the mode-coupling functional, a quadratic polynomial that can reproduce all vahd values 
for the exponent parameter A is used |pi||: 



m{t) = v,(Pit) + V2 {(p{t)f . (3b) 

For A = 0, ideal liquid-glass transitions occur on a line in the V1-V2 plane of coupling 
constants. One can use A to parameterize this line of critical coupling constants: 



v1 = (2A - l)/A^ = VA^ 1/2 < A < 1 . (4) 

Thus this model is specified by two control parameters (fi,f2), by two frequencies 
quantifying the transient dynamics, and one rate A for the activated transport processes. 
The model has many non-generic features, and therefore one cannot expect it to describe a 
measured spectrum. In the present paper, the correlator 0(t) is introduced to mimic in an 
overall fashion the combined effect of all structure fluctuations in producing the bifurcation 
point and the exponent parameter A of the system. 

The dynamics of some probing variable A coupling to density fluctuations shall be de- 
scribed by a second correlator, to be denoted 4>\(t). It obeys an equation analogous to 
Eq. (P: 

<pAi^) = ... ^.,,r~^ . . .1 • (5a) 



Again the microscopic dynamics is quantifled by two frequencies referred to as microscopic 
parameters (f2^,z/^). The activated relaxation processes are described by A^. The mode- 
coupling functional shall be specifled by a coupling to (pit) quantifled by a single coupling 
constant v\: 

m\{t) = v\mm) ■ (5b) 

It is a peculiarity of this model, that the dynamics of the probing variable A is influenced 
by 0(t) but not vice versa. Thus the position of the transition is not modifled by the 
introduction of the second correlator nor is the value of A. The model was motivated by 



Sjogren for the description of tagged-particle motion in a glassy environment, and it will 
be used here in the same context for the interpretation of the neutron-scattering data. The 
MCT for the reorient at ional dynamics of a non-spherical probe molecule suggests the same 



schematic model for the dipole and quadrupole relaxation |^ ; an observation that motivates 
the application of the model for the description of the dielectric-loss and depolarized-light- 
scattering spectra, respectively. For the incoherent-neutron-scattering cross-section the flt 
will be done using the model parameters for 0^(t) different for different wave vectors. For 
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the index A the abbreviations Is, de, and ns for hght scattering, dielectric loss, and neutron 
scattering, respectively, will be used. The specified two-component schematic model has 
been used earlier for data interpretation with the restriction to A = = 0. Depolarized- 
light-scattering spectra within the full GHz band have been described for glycerol for all 



temperatures above Tg [^, and for ortho-terphenyl for T > |]2^. Ruffle et al. were 
the first to simultaneously describe glassy spectra for several probing variables A. Within 
the /3-relaxation regime, they fitted coherent-neutron-scattering spectra for several wave 
vectors and also the longitudinal elastic modulus for Nao.5Lio.5PO3. 

The single coupling constant v% determines all features of the structural-relaxation part 
of the second correlator. Thus, the a-peak strengths, widths, and positions are correlated. 
These correlations follow the same pattern as found and explained for the a peaks of the 



hard-sphere system |]T3|,|T^. Nevertheless, it is not obvious from the beginning, and thus 
truly remarkable, that such a simple model will be sufficient not only to explain the trends 
found in the data, but even to reproduce structural relaxation for PC quantitatively. 
Equation (^) is equivalent to 

4>{t) + (A + z/)0(t) + {n^ + A- u)^{t) 

f m{t - f) U(t') + A ■ (pit')] dt' = , (6) 

Jo L J 

to be solved with the initial condition </)(t = 0) = 1, </)(t = 0) = — A. This equation, together 
with Eq. (|3ED, is solved numerically with a similar algorithm as used in the preceding work 
for the case A = 0. Equations (^) are treated in the same manner, but (f){t) has to be used 
as input for Eq. (|5b|) . From the result for a Laplace-transformation yields (p^iz). 

The fluctuation-dissipation theorem then determines the dynamical susceptibility XAiz) of 
variable A: 

Xa{z)/xa = z(t)\{z) + 1 . (7) 

Here, xa oc {A^) is the thermodynamic susceptibility. In particular, the imaginary part 
of Eq. d^) determines the normalized susceptibility spectrum, x"Ai'^)IXA = ^^a"!"^)? 



quantity of main interest in the following. In our data analysis, xa enters as an additional 
fit parameter, which we treat, for the sake of simplicity, as a temperature-independent 
normalization constant. 



III. DATA ANALYSIS 



A. Fits to the data 



The result of our fits to the measured PC spectra are shown by the full lines in Figs. || and 
^. Since one cannot expect the schematic model to provide a description of the microscopic 
band, the fits have been restricted to frequencies below 500 GHz for the light-scattering and 
neutron-scattering spectra. The fit range for the dielectric spectra could be extended up 
to 1 THz. For the neutron-scattering data, a set of spectra for 3 representative g-vectors 
out of 10 analyzed is shown. The analyzed g-range is 0.5 A < g < 1.4 A ; outside this 
range, experimentally accessible frequency windows become too small to gain meaningful 
information for MCT parameters. In Ref. [|l|, light-scattering spectra above T = 250 K have 
been published, but show apparent violation of a scaling. We were able to fit these curves 
with the same quality as the ones shown by assuming a slightly varying static susceptibility 
Xis, which has the effect of shifting curves up and down in the log-log plot. These curves 
were omitted in Fig. |T] to avoid overcrowding. 

All model parameters should be used as temperature-dependent fit parameters in our 
analysis. Within the studied temperature interval, there are no structural anomalies reported 
for PC. Thus, the fits are done with the constraint that the parameters drift smoothly and 
monotonously. In the following part of this section, the parameters used for the theoretical 
curves in Figs. |l] and shall be discussed. 

One experiences a considerable fiexibility in choosing the path (wi(T), t>2(T)) followed 
by the coupling constants in the f i-f 2-parameter plane for the interpretation of the data as 



emphasized earlier ||2^. To arrive at an overview of the possibilities for fitting the many 
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spectra, we started with a first step, where the path was varied but biased to some smooth 
curve. Applying the general theory to Eq. one derives the formula for the above- 
mentioned separation parameter a, 

a = (1 - n [{v, - vi)r + {v^ - v'i)r^] . (s) 

In our first step of the analysis, we also force the fi, V2 to obey the asymptotic linear 
(Tc — T) dependence of a cited above. In the second step, this latter restriction is eliminated 
and a free fit is started by examining small corrections to the result of the first step. The 
thus obtained results also account for an inevitable uncertainty in the determination of the 
experimental temperatures. The fit yields Tc ~ 180 K, and A ~ 0.75, corresponding to 



a ^ 0.30 and b ^ 0.56. The value for A is between the values reported in Refs. [|1| and |T0 
and falls within the error bars of both. The linear interpolation of the found a versus T 
values gives a = C(Tc — T)/Tc with C ~ 0.069. The found distribution of (^1,^2) points 
is shown in the upper part of Fig. |^. Upon lowering T, both vi and V2 increase, which is 
consistent with the physical reasoning of the system's mode-coupling coefficients becoming 
larger at lower temperatures. The lower diagram in Fig. ^ demonstrates that the asymptotic 
formula for a is well obeyed for 150 K < T < 285 K. It should be stressed that the glass- 
transition line is just crossed by a regular drift, i. e. there is no accumulation of (fi, V2) points 
close to it. This demonstrates how the critical phenomena predicted by the MCT originate 
from the mathematical structure of its equations of motion. In particular, the schematic 
model illustrates that within MCT no subtle g-interferences or hydrodynamic phenomena 
are responsible for the glass-transition dynamics. 

The fitted mode-coupling coefficients v%{T) for the light-scattering and dielectric data, 
and the corresponding coefficients v^^{q, T) for the neutron-scattering experiment are shown 
in Fig. ^. Again, we find monotonically increasing couplings with decreasing temperature. 
The coupling coefficients v^^{q) describing the incoherent-neutron-scattering data are de- 
creasing with increasing q. This is equivalent to the plateau values f^''^ decreasing with 
increasing q, which agrees qualitatively with the findings for incoherent-neutron-scattering 
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results discussed within the microscopic MCT [0 . 

The parameters fi^, which specify the transient dynamics of (f)\ are shown in Fig. |^. 
The resuhs from the neutron-scattering analysis reflect the behavior fi^g(g,T) oc q ■ \/T to 
a good approximation, which is in agreement with the result of the microscopic theory. But 
drawing more conclusions from the microscopic parameters would be over-interpreting the 
model. They are shown here mainly to demonstrate that there are no abnormal variations 
occurring. We find much larger uncertainties for the microscopic fit parameters f2, z/, f]^, 
z/^, than for those parameters fi, f2, and v% ruling the structural-relaxation part of the 
spectra. In particular, it was possible to use for the parameters which specify the transient 
of the first correlator temperature independent values = 1 THz and i/ = THz. 

The hopping coefficient A in Eq. (|3a|) determines the position of the susceptibility min- 
imum below Tc. This minimum cannot be seen in the light-scattering data, thus the chosen 
values are not unambiguously determined. The light-scattering spectra in the upper panel of 
Fig. |I] are fitted with the hopping parameter for the second correlator ignored: Aj*g = 0. 
The fits to the dielectric-loss spectra in the lower panel of Fig. [1| are done with a non- 
vanishing Ajg. For the whole temperature range investigated, A(T) can be assumed to 
follow an Arrhenius law, A(T) oc exp(— E'^/T), which would be expected for thermally acti- 
vated hopping over barriers. Fig. ^ shows the values used for the fit. Although A increases 
by an order of magnitude, the calculated curves for temperatures higher than 190 K show 
no infiuence from hopping effects on the spectra. This is demonstrated in Fig. |^. The 
irrelevance of the increasing hopping coefficients A^ for temperatures increasing above 
can be understood on the basis of a discussion of the asymptotic formulas [^]. It is the 
reason, why the idealized theory can be used for data analysis for T sufficiently larger than 
Tf.. In the analyzed neutron-scattering experiment, the dynamical window and the studied 
temperature intervals are too small to investigate hopping effects, and therefore the curves 
in Fig. are calculated with A^^ = 0. 

Above Tc, the spectra including hopping show deviations from the idealized ones only 
for small T — T^. Below Tc, the crossover to the white-noise spectrum is suppressed, and 

12 



a minimum occurs as hopping starts to be the dominant relaxation effect. Because of 
the insensitivity of the main body of the analyzed data to choices of A, the activation 
energy cannot be determined very precisely from the fit; the upper straight line in Fig. |^ 
corresponds to Ea = 811 K. This value is in reasonable agreement with the one found 
in an earlier asymptotic analysis []T[. Dielectric- loss spectra show hopping- induced minima 
at higher frequencies than the light-scattering spectra, and this we have accounted for by 
introducing a second hopping parameter A^g there. In a similar way, (fi^g/f2)^Ajg(T) follows 
an Arrhenius law and has no influence on the spectra above T^, this second hopping term 
has already been included in the comparison studied in Fig. |^. Here, the activation energy 
is of the order of 2000 K, which makes the result more striking, since (f2^g/fi)^A^g is allowed 
to vary over three orders of magnitude. In both cases, activation energies as well as the 
prefactors are of reasonable magnitude. It should be stressed that, although the treatment 
of hopping by a frequency-independent A is rather crude, the resulting frequency range in 
which the schematic model gives a good fit to experimental data, is enlarged by about one 
decade for T < relative to the fit interval which can be treated by the idealized-MCT 
model. 

In the measurements of the dielectric functions, information on both the imaginary and 
the real part of eiu) = e'{uj) + ie"{uj) have been obtained The fit to the e" data 
shown above was performed using e"{uj) = Airx'^^iuj) = 4:TcXdc^4>de"i^)^ ^^us obtaining the 
proportionahty factor Eq = 4:nxde as a by-product. Then, the real part is given by s'{uj) —i = 
Eq ■ {1 + u;0''^g(u;)). The new parameter i has to be determined by shifting the curves, and 
it can differ from ^oo = 1 in both directions: The liquid exhibits microscopic oscillations, 
which contribute to e'{uj) as some shift Ae^j^j, with respect to £00 = 1 for the structural 
part of the response function. The schematic model uses a single damped oscillator, giving 
some A£:^jj.j. , which may be either too small or too large. Depending on the temperature, 
we find values of e = £00 + (Ae^cr. ~ ^^micr.) between 3 and —1, which are of reasonable 
magnitude. Figure |^ shows the result of testing our fit against the accordingly shifted 
real part of the measured dielectric function. It is clear from the theory that the real and 
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imaginary parts of the calculated curves are connected by Kramers-Kronig relations. But 
for the experiment, both quantities have to be regarded as almost independent data sets, 
since the measurements are restricted to a finite frequency range. Thus, figure § provides 
more than just a different view on the fit shown in Fig. |1|, and it is an important point 
that the real-part data can be fitted with the schematic model as well, introducing only 
one additional fit parameter i. In the minimum region of the spectra, we find this to be 
confirmed, and for higher T, the a-relaxation step can be described by the schematic model, 
too. The discrepancies for the a peak in the glass are the analogue to what can be seen 
in the e" fit. Similar observations hold for the high-frequency dynamics, where one has to 
notice in addition, that experimental error bars are relatively large for frequencies above 
300 GHz. A slightly better fit of the e' data could have been achieved by allowing the static 
susceptibility Xde to vary with temperature. This possibility is not examined here, since the 
shift is only small, and since we do not want to introduce assumptions on the T-dependence 
of the static quantity ^o- 



B. Summary of the Data Analysis 

Glass-forming liquids exhibit temperature-sensitive spectra for frequencies well below 
the band of microscopic excitations. These precursors of the glass transition are referred 
to as structural-relaxation spectra. The full lines in Fig. || and demonstrate, that the 
evolution of structural relaxation of PC, including the crossover to the microscopic regime, 
is described well by a schematic MCT model. The description holds for all spectra obtained 
by the depolarized-light-scattering spectrometer; in this case it deals with the three-decade 
dynamical window between 0.3 and 500 GHz, and it accounts for the change of the spectral 
intensity by a factor of 10^ if the temperature is shifted between the glass transition Tg and 
30 K above the melting temperature T^. It accounts for the measured a-peak-maximum shift 
by a factor of 10 if T is changed by 30 K. A similar statement holds for the description of the 
dielectric-loss spectra, where the a-peak shift from 40 GHz down to 0.02 GHz is described. 
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This shift is caused by a temperature decrease from 293 K = + 75 K to 243 K. 

Between the a peak and the vibrational excitation peak near 1 THz, the susceptibihty 
spectra in Figs. |1| and exhibit a minimum at some frequency cUmin- It shifts to smaller 
frequencies as the temperature is lowered, but less than the a-peak position. Its intensity 
Xmin = x"(^min) cxcccds the white-uoisc spectrum one would expect for the dynamics of 
normal liquids by more than two orders of magnitude. Such white noise would yield suscep- 
tibility spectra varying linearly with frequency, Xwh.n.('^) '^i is indicated by the dashed 
lines in Fig. |l|. These anomalous minima are also treated properly by the model. 

Neutron-scattering data are available for a series of wave vectors g, and hence the dynam- 
ics is probed on various length scales. The q dependence is in the schematic model described 
by that of the coupling coefficient fns('?)- The data description in Fig. ||is possible using a q 
dependence in qualitative agreement with the results expected from the microscopic theory 
of simple systems. 

It appears nontrivial that the used schematic model can deal with the mentioned spectra 
of PC. The success of the fits indicates that the studied glassy dynamics is rather insensitive 
to microscopic details of the systems. Apparently the evolution of glassy dynamics within 
the GHz window refiects, above all, only quite general features of the nonlinear-interaction 
effects, which can also be modelled by simple truncations of the full microscopic theory. 
These conclusions require some reservation. The explanation of the PC data by the used 
model is based on the choice of the model parameters, in particular on the choice of the drift 
of all parameters with changes of temperature, which is documented in Fig. Only a full 
microscopic theory can show, whether or not the chosen parameters are in accord with the 
fundamental microscopic laws. 

Furthermore, it has to be emphasized that the studied model cannot reproduce the 
spectra for frequencies below 1 GHz if the temperature is below the critical value T^,. Such 
spectra can be measured accurately using dielectric-loss spectroscopy, and the lower panel 
of Fig. |I| exhibits some of this data for T = 173 K and T = 183 K. The lack of success 
of our work in handling these spectra is clearly connected with the improper treatment 
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of hopping processes. It remains unclear at present, whether this is due to the stochastic 
approximation, Ng{z) = iAg, or due to restricting ourselves to a one-component schematic 
model, or whether the whole extension of MCT to a theory including hopping transport is 
inadequate. 



IV. SOME ASYMPTOTIC FORMULAS 

Let us list some of the asymptotic results for the studied MCT model, which will be 
needed below in Sec. 0. These results are obtained by straight-forward specialization of the 
general formulas discussed in Ref. [^. We will focus on the /3-relaxation regime for T > Tc, 
with hopping effects neglected. A comprehensive discussion of the asymptotic results can 



be found in Ref. 13 



From the full MCT equations ([1|), a leading-order expansion in y |o"| gives rise to the 
asymptotic predictions for the intermediate-time window of the (3 relaxation. A central 
result is the factorization theorem, — fq = hqG{t), where the so-called j3 correlator 

G(t) is independent of q. This result still holds, in the generic case, for the tagged-particle 
density-fluctuation correlator or the correlator dealing with light scattering or dielectric 
response: = f^'^ + h\G{t), with the same G{t) as above. The Fourier-cosine transform 

of G{t) is called the /5 spectrum G"{uj). One gets for the normalized susceptibility spectra 

x:{oo) = u;<j>:{u;) = K^\uj) , (9) 

where = ojG"{ijj) is called the /5-susceptibility spectrum. Here, the index x denotes 

either the wave-vector modulus g, or x = (s, A). The function G depends on t/to, cr, and 5 
only: it is uniquely determined by the exponent parameter A as the solution of the equation 

a-5t + \{G{t)f = ^ f G{t - t')G{t') dt' , (10) 
at Jo 

to be solved with the initial condition G{t — > 0) = (t/to)""". The so-called hopping parameter 
6 has to be calculated from A^, and for the studied model it reads 
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5 = Ar7(i-r). (11) 

In this context, the numbers only enter as corrections to scahng. 

The plateau values and the critical amplitudes can be calculated from the mode- 
coupling functionals. In the case of the schematic model studied, the values for the first 
correlator are given by A: 

r = i-A, h = {i-n. (12) 

The relation between the exponent parameter A and the a-peak strength is one of the 
non-generic features of that model. For the second correlator, the plateau value and critical 
amplitude read 

1 I - f 

f/ = 1 _ h'. = —\ . (13) 

Changing v%, the a-peak strength f^'^ can be varied. Again, these equations establish a 
non-generic relation between the f^'^ and the h\. In our fits to the neutron-scattering data, 
a q dependence of flf and k^^ can arise only through a g-dependence of the v^^. 

From Eq. (0) one identifies for the case 6 = the time scale for the /3-relaxation: 
ta = to\a\^^^'^"-. Going over to rescaled times, i = t/t^, and rescaled frequencies, Cj = tota-, 
one gets from Eq. @ the scaling law for the /^-susceptibility spectra 

x"{uj) = KcaXii^) , (14) 

where Co- = \f\o'\. The master spectrum x is cr- independent. It is fixed through the exponent 
parameter A, and thus through the static structure alone. For large rescaled frequencies, 
ci) 1, one obtains the critical-power-law spectrum. This extends to all frequencies as 
a ^ 0: 

x'iu) = h,-sm{7ra/2)Til-a){ujtor, T = T, . (15) 

For small rescaled frequencies, one gets the von Schweidler-law for a < 0, x(cj ^ 1) oc 
and thus x exhibits a minimum at some frequency cjmin with Xmin = x(^min)- Due to 
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the scaling law, Eq. (0), the variation of the spectral minima with temperature is, in the 
asymptotic region, given by 

^min ^min/^(T5 Xmin Xmin ' CT < . (-^^) 

The point ((^min, Xmin) is completely fixed by A, and for A = 0.75 one gets: cDmin = 1-733, 

Xmin ~ 1.221. 

On the glass side, a > 0, the idealized theory yields for the (3 correlator for large rescaled 
times a constant, G(t ^ 1) = 1/^/1 — \. Thus the signature of the MCT-fold bifurcation 
are y/Tc — T anomalies of the nonergodicity parameters fx = (pxit oo): 

f,{T) = f: + K^a/{l-X), T<n. (17) 

If the correlators deal with density fiuctuations or tagged-particle densities, the quantity fx 
is the Debye- Waller factor or Lamb-MoBbauer factor, respectively. For a < 0, corresponding 
to T > Tc, the long-time limits of the correlators vanish, as is the case for T < Tc but 6^0. 
But if a and 6 are sufficiently small, the correlators still exhibit plateaus for times exceeding 
the transient scale to before the decay towards zero sets in. The heights of these plateaus 
are then given by fx for T < Tc, and by for T > Tc, then called effective nonergodicity 
parameters. The decay from the plateau is the a process, and thus the strength of the a 
peak in the susceptibility spectra is given by fx- This also corresponds to the height of the 
relaxation step exhibited by the real part of the susceptibility, when the frequency is shifted 
through the a-peak window. 

The preceding Eqs. (P-|T7|) establish universality features of MCT. They provide the basis 
of a general explanation of the glassy MCT dynamics by means of features of the spectra 
not depending on the specific microscopic properties of a given system. 

V. SCALING LAW ANALYSIS 

In this section it shall be studied how well the above calculated MCT solutions can be 
described by the MCT-/3-relaxation-scaling laws summarized in the preceding section. It has 
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been demonstrated earlier [T^JTJ], that the range of vahdity of these equations can be ana- 
lyzed by evaluating the next-to-leading-order corrections. Here, we will study the combined 
effect of all corrections due to structural relaxation as well as due to vibrational-transient- 
dynamics effects. Only the solutions referring to the parameter sets used in Figs. |I] and |^ 
will be discussed. Thus the following analysis refers to control parameters and dynamical 
windows representative for state-of-the-art experimental studies of the evolution of glassy 
dynamics. 



A. The critical decay 

Solving the equations of motion for T = Tc, A = 0, and A* = for times up to 10^^ ps, 
the critical power law, Eq. (H), was identified. The common time scale was determined to 
to = 0.035 ps. The leading-order result ^^(t) = (t/to)~", where (l>xit) = {(pxit) — fx)/f^x, is 
shown in the double-logarithmic representation of Fig. ^ by straight dash-dotted lines with 
slope —a. For the two temperatures closest to Tc, the full lines in this diagram exhibit the 
solutions (f)x{t)- Dashed lines demonstrate the corresponding (3 correlators G{t), determined 
from Eqs. (H,|lO|,|lT]) . The approach of the first correlator towards the plateau is well 
described by the scaling law for T = 180 K and 190 K. For T = 180 K the critical power law is 
exhibited within a 1.5-decade time window for times exceeding with tc/^o ~ 300, while for 
t < tc the vibrational transient dynamics masks the structural relaxation. In this case, the 
validity of the critical power law for larger times is restricted by the onset of hopping effects. 
Hopping plays no significant role for the (3 relaxation of T = 190 K (compare Fig. |^. But 
there, the deviations of G{t) from the short-time limit (t/to) " set in already for t < t^. Thus 
this power law cannot be identified anymore for distance parameters e = {Tc — T)/Tc with 
l^l ^ 0.06. This scenario is in semiquantitative agreement with the one discussed in Ref. |]2^ 
for the density correlators of a hard-sphere system. Let us reiterate that the correlator (f){t) 
drives the glass transition for the studied model, but that it is not the quantity measured. 

The two lower sets of curves in Fig. |^ show that the decrease of 0fs(t) and 0de(^) towards 
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their plateaus /ig'*^ and f^^ respectively is described qualitatively by the dashed lines, i. e. 
by the scaling laws. However, there are remarkable quantitative deviations between the 
solutions (l>%(t) and their asymptotic form Git). These appear as if the amplitude experiences 
some offset. The reason is that the transient dynamics influences the correlators (f)%{t) also 
for times which exceed tc by up to two orders of magnitude. This means that the dynamics 
of the two probing variables is strongly influenced by oscillations within that window, where 
the driving correlator exhibits the law. Therefore the power law = (t/to)~" 

cannot be identified accurately in the curves shown for A = Is and A = de. This is also 
demonstrated by the straight dash-dotted line in the upper panel of Fig. |l|, which represents 



the asymptotic low-frequency-susceptibihty spectrum at the critical point, Eq. (pISl). 

Within the 1.5-decade window, where the log0-versus-log(t/to) curve for 180 K in Fig. |^ 
demonstrates the critical-decay asymptote, the graphs of log0fg and log^^g versus log(t/to) 
for 180 K and 183 K, respectively, also appear as nearly straight lines, so that they can 
be described very well in this window by some effective power law. One thus expects an 
effective power-law spectrum which is described by Eq. (|15|), but with a and Ha replaced 
by some a^^ and /i^ respectively. This phenomenon also was observed for the susceptibility 



spectra of the hard-sphere system For the light-scattering result, one infers a < a 

and hf^ < hf^. The dotted line in the upper panel of Fig. |l] corroborates this conclusion. It 
exhibits the solution for the model evaluated for T = with hopping effects ignored. This 
line can be fitted well between lO^^THz and lO^^THz by an effective power-law following 
Eq. (|15D with a^^/a ^ 0.92 and h^^/h ^ 0.7. The crossover from this effective power law to 
the asymptotic critical law, Eq. (|1^), occurs only at frequencies around 1 MHz. 

B. The non-ergodicity-parameter anomaly 



Figure |T^ shows effective nonergodicity parameters of the three correlators underlying 
the curves in Fig. [1|, which were determined from the plateau heights of the 0x(^)-versus- 
logt diagrams. The crosses in the lower panel show the values deduced in Ref. M from 
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the step size of the measured real part of the dielectric function, divided by the value of Eq 
assumed in the fit to the susceptibility spectra. Figs. |I] and |^ demonstrate that the present 
model describes the dielectric function of PC reasonably, and so it is not surprising that 
the calculated values (dots) reproduce the measured ones (crosses) reasonably well. The 
discrepancies between dots and crosses are anticipated to be mainly due to difficulties in 
determining the step size accurately in the experiment, where one carefully has to eliminate 
contributions from the f3 relaxation. 

Full lines in Fig. |10| exhibit the asymptotic laws, i. e. the values fx from Eq. ([17D for 
T < Tc and the constant for T > Tc. Figure |lO|(a) demonstrates that the 60% variation 
of the effective non-ergodicity parameter / of the first correlator is described well by the 
asymptotic formula. This holds for temperatures down to Tg. On the other hand, the results 
for the light scattering and for the dielectric response do not exhibit the asymptotic behavior; 



there is no evidence for the a/Tc — T anomaly at all to be noticed in the data. There are 
two reasons for this finding. The obvious one refiects the large size of f^'^, i. e. it results 
from the observation that the a peaks of the susceptibility dominate over the remaining 
susceptibility spectrum (compare Fig. [1|). Equation (p!7| ) for probing- variable correlator is 
equivalent to 

(1 - fl) = (1 - r/) - h\^a/{l-\) . (18a) 
Since (1 — /^) and h\ are positive and (1 — f^/^) is less than 0.1 for the two correlators 



discussed, the whole a/Tc — T effect is below 10%. Therefore it is difficult to separate the 



a/Tc — T anomaly from the scatter of the data. The less obvious reason results from the 
smooth but appreciable temperature drift found for the coupling coefficient v\ (compare 
Fig. This coupling determines the non-ergodicity parameters of the second correlator of 
the schematic model in terms of the parameter /: (1 — /^) = l/{v\f). The square-root 
singularity is due to that in /, and expanding (1//) one reproduces Eq. ( 18a ), but with 
effective terms 

(l-/l'T'^ = i?^(l-/r), h^^ = RAh\, 
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Ra{T) = vr/v\{T) . (18b) 

Replacing the renormalization coefficient Ra{T) by its value at the critical point, = 1, 
one reproduces the leading-order result, Eq. (|1^. However, within the temperature interval 
considered, the smooth drift of (1 — f^'^Y^ overwhelms the small variation of the y^Tc — T 
term. This is demonstrated in Figs. ^OKb,c) by the dashed lines. The numerically found 
circles are well described by this line. One concludes that the drifting coupling coefficient 
v'\{T) is responsible for the deviations from the leading-order asymptotics. Unfortunately, 
the v%{T) are not available directly from experiments. 



C. Scaling of the /3-relaxation minima 



Figure |ri| shows the susceptibility master spectrum x for A = 0.75 and 5 = as dashed 
curves. The upper set of solid lines in this figure are the spectra of the ffist correlator 
rescaled, according to Eq. (p!^), as u(f)"{uj)/{J\a\h). Asymptotic validity of scaling is demon- 



strated: the window of rescaled frequencies uj = ujt„-, for which the rescaled spectra are close 
to the master spectrum x, expands with decreasing (T — T^. Convincing agreement be- 
tween X and the 180 K result can be found as long as hopping effects are ignored. For higher 
temperatures, where \e\ = \T — T^jTf. > 0.06, strong deviations are found. The T = 210 K 
spectrum, for which \e\ = 0.17, does not even show a minimum. The demonstrated devi- 



ations from the scaling laws are similar to what was explained in Ref. [T^ for the MCT 
solutions for the hard-sphere system. 

Preasymptotic-correction effects for the variables discussed for PC in Figs. |I| and differ 



from those for the auxiliary correlator (j). This is demonstrated in the lower part of Fig. |Tl 
for e = —0.17. Deviations of the rescaled spectra x'a{^) /(V^^a) from the master spectrum 
x{ujt„) are larger for the dielectric loss than for the light scattering, and the latter are larger 
than those for the neutron-scattering results. While the predicted probe independence of 
XA{^)/hA holds rather well for u < ujj^i^, deviations from the factorization theorem are 



observed mainly for higher frequencies. As discussed above in connection with Eq. (|18aD 
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the large size of f^"^ leaves only a 10% decay of the correlator from the initial value unity 
to the plateau, and this decay is influenced by vibrational motion. This leads to the strong 
disturbances of the susceptibility spectra for uj > cUmin- For the neutron-scattering data for 
intermediate wave vectors, this problem is not so severe, since the critical Lamb-MoBbauer 
factor decreases with increasing q. Therefore the shape of the susceptibility minimum, 
which is exhibited by the two neutron-scattering results shown in Fig. 0, is closer to the 
one of the master spectrum. 

The g-dependence of the critical amplitude for the incoherent-neutron-scattering 
spectra has been measured as a byproduct of the test of the factorization theorem, 
(pq'iyj) oc h'^^'x" {oj) / oj . A linear law, /i^ oc g, was found within the studied wave- vector 
interval [1^. Such a strictly linear law is not compatible with the microscopic MCT, which 



predicts that the /ig-versus-g graph exhibits a broad asymmetric peak near the position g^ax 
of the first sharp diffraction peak of the structure factor. For small g, the critical amplitude 
increases regularly as hq + 0{q'^) |T^, and thus hq exhibits an inflection point for 

some g < gmax- Whether the found linear g-dependence of the experimental values is due 
to multiple-scattering effects, is not clear |T^. Our schematic-model fit, however, suggests 



that, even if can be approximated by a linear law, the g-dependence is not strictly linear 
but rather given by some intermediate crossover around the inflection point. Equation ([13D 
relates /i^ to the inverse of f^; thus a strictly linear law for /i^ would imply oc 1/g. Such 
result is added as a dash-dotted line in Fig. ||b, and it shows that this is not consistent with 
our data analysis. An ad-hoc expression, reflecting the crossover from the small-g asymptote 
through the inflection point is oc g^/ [1 -|- (g/g*)]. The resulting expression for is added 
in Fig. ^3 as dashed lines for two temperatures, and it provides a reasonable interpolation 
of the found fit parameters. One would need measurements for wave vectors of the order 
of 0.2 A ^ and less, in order to test for the small-g^ behavior predicted by the microscopic 
MCT. 

There is a most bothersome preasymptotic-correction effect which can be seen in the 



lower part of Fig. [IT]: the positions i^anata- of the susceptibility minima are not identical; 

23 



and they are all larger than the asymptotic value cjmin which is shown by a diamond. Since 
the bands of microscopic excitations of the correlator spectra (j)^'^{uj) are located at much 
lower frequencies than that of the spectrum (f)"{uj), the spectrum of the test variables cross 
over too quickly to the transient to be able to develop the universal relaxation pattern for 
u > cjmin. As a result, tUmin gets an offset to larger frequencies. Figure O exhibits this 



result as a rectification diagram. The asymptotic result is shown as a full straight line: 
'^min = (^^minAo)^1c^| = C{T - with the constant C = (cl^mm/to)^" ■ C'/Tc ^ 0.004. The 
positions of the observed minima can still be interpolated reasonably by straight lines, shown 
in dashed. However, the slopes of the dashed lines differ from those of the asymptotic line. 
In a clear violation of the asymptotic factorization theorem, the lines for different probing 
variables are different. The linear interpolations lead to intersections with the abscissa, 
which differ somewhat from the correct value of T^. For the neutron-scattering data, this 



interpolation has been omitted in Fig. |T2|, since the error bars obtained by g-averaging do 
not allow for a well-determined estimate here. 

There is also a strong temperature-dependent offset of the amplitude scale relative to 



the scaling-law prediction. This is demonstrated in Fig. O for the light-scattering spectra. 



All four rescaled minima Xmin/v^^is are far below Xmin, which is indicated by a diamond. 
Moreover, with decreasing (T — Tc), the discrepancy between rescaled curves and expected 
asymptote does not decrease, rather it increases. Such behavior is not anticipated from the 
leading-order corrections to the scaling laws ||T3| , p!^ , but it can be explained as a higher- 
order effects because of the important role played by the temperature dependence of the 
coupling coefficient f^. This drift can be eliminated by introducing an effective amplitude 
hf^, as discussed above in connection with Eq. ( p.8b|) . The result is given by the upper set 
of curves in Fig. [1^. Indeed, the discrepancies between asymptotics and rescaled curves are 
reduced, and they now decrease with decreasing (T — Tc). But even for T = 190 K, i. e. for 
l^l = 0.06, there is a considerable offset of the minimum intensity from the scaling result. 
The 180 K curve demonstrates the approach towards the asymptotic limit, would hopping 
be absent. There still is a clear deviation between rescaled spectrum and scaling-law result, 
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which increases with increasing uj for uj > cUmin- But the sign and size of this effect are 
similar to what was found for the hard-sphere system for wave vectors yielding a plateau 
as large as ff^ |T^. 



D. Summary of the Scaling-Law Analysis 

It is, of course, more satisfactory to interpret data for glassy dynamics with the set of 
universal formulas provided by MCT for the asymptotic dynamics near a glass-transition 
singularity than to explain experimental findings within schematic models. The more prob- 
ing variables A are taken into account, the more convincing such an analysis is, since the 
universal results also imply connections between spectra measured for different A. The pre- 
ceding work on PC [|l],^-^|10| exemplifies these statements. However, the data are influenced 
by preasymptotic effects, and one cannot judge the relevance of these correction effects, if 
one does not know the underlying microscopic MCT equations. Forcing data into the uni- 
versal formulas can thus lead to self-contradicting results, as the preceding subsections have 
demonstrated. While the spectral shapes are rather robust and the rectification diagram 
for the scales appears correct and leads to a reasonable estimation of Tc, as shown by the 



dashed lines in Fig. |12|, the prefactors for the asymptotic formulas extracted from the data 
can be quite wrong. This error cannot be noticed if one studies a single probing variable A 
only, but it appears as a violation of the factorization theorem if one compares spectra for 



different A. One concludes that the problems with the analysis discussed in Ref. |]T0[ are 
neither due to inadequate application of MCT results nor due to failures of MCT. Rather 
they reflect the properties of MCT; more precisely, they exemplify the limitations for the 
application of asymptotic laws. 



A general rule [0,0 for the test of the /3-relaxation-scaling law is corroborated by the 
present analysis: if the nonergodicity parameter is large, i. e. if the a-peak strength 
is large compared to the strength of the microscopic-excitation peak of the susceptibility 
spectrum, the preasymptotic corrections are very important. This is especially true for the 
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discussed light-scattering and dielectric-loss spectra. Neutron-scattering spectroscopy has 
the advantage that can be shifted by changing the wave vector q. Therefore, we found 
the scaling-law analysis to work best for the neutron-scattering data of Ref . ||10| . It would 
be very informative to corroborate this finding by a measurement of the expected — T 
anomaly of the Debye- Waller or Lamb-Mofibauer factor. 

VI. THE DIELECTRIC MODULUS 

Since the memory kernel m'^{t) in the mode-coupling approach is expressed as a polyno- 
mial of the density correlators, Eq. (p^) , this quantity shows the same asymptotic scenario as 
the correlators themselves. From the factorization theorem for the correlators one concludes 
in leading order for the kernels: = f^^j^ + ^G(t). Eqs. ( ^ , p!2|jl3|) determine the 

plateau flfj^ and the critical amplitude h^^ ^ for the memory kernel of probing variable A: 

Jm,A — ^AJ J a 5 

hM,A = vl irh\ + fTh) = v\h . (19) 

While in the above discussion f^/" was found to be larger than 90%, such that the square- 
root singularity is suppressed to a below-10% effect, the situation for the memory kernel is 
different. The coupling coefficient v\ now plays the role of a normalization constant. If one 
introduces the normalized memory kernel in analogy to the normalized correlators (pxit), 
^a(^) = ^Ai.'t)/'^A^ such that m\{t ^ 0) = 1, one gets for the normalized plateau: 

fi:U = rr/ ■ (20) 

In cases where /^''^ is close to unity, like in our analysis of the dielectric-loss and light- 
scattering data, one can approximate m^(t) — + hG{t). This equals the asymptotic 
expression for the first correlator. Thus, we can expect (3 scaling for the memory kernel of 
the probe- variable A to work equally well as for the first correlator and thus better than for 
the corresponding probe-variable correlator. Let us examine this in detail for the memory 
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kernel of underlying the fit to the dielectric-susceptibility spectra. For the light-scattering 
data, qualitatively the the same picture arises. 

In the upper part of Fig. |14|, three spectra of the memory kernels for the dielectric 
function, rescaled to —ujm'^^{uj) / ^Jlo], are plotted as solid lines for three temperatures above 
Tc- The asymptotic /5-susceptibility spectrum for A = 0.75 is again shown as a dashed line. 



While the picture shows some similarity to the situation found in the upper part of Fig. ^ 



for the light-scattering susceptibilities, the reason for the deviations from the scaling law are 



different. This can be inferred from the lower part of Fig. where the same scaling is shown 
for the normalized memory functions. Here, the solid lines represent —ojm'^^iuj) / -y/jo^/i. One 
notices, that the a-peak strength is remarkably smaller than in the dielectric susceptibility, 
and comparable to that for the first correlator of the model. Similarly, we find the standard 
scenario for the approach of the rescaled spectra to the master curve. The deviations from 



the asymptotics are qualitatively the same as exhibited in the upper part of Fig. 11 for 



uj(j){uj). Thus one concludes: the deviations from scaling seen in the upper part of Fig. |n 
are mainly due to the T-dependent normalization and not, as in the case discussed in 
connection with Fig. |T3|, due to microscopic crossover effects. 

The question of normalization becomes even clearer for the effective non-ergodicity pa- 
rameters. Figure |l^(a) shows the unnormalized values /X/je open circles. The full 



line exhibits the asymptotic prediction /X/^e — fuAe + ^M,de\/^/(l ~ ^) f'^^ T < Tc and 
flide — fhide ^ — ^c- Again, the drifting coupling coefficient is responsible for 



masking the predicted square-root law. But, unlike in Fig. [1^, this is only true for the 
unnormalized quantity. The normalized function flj^g, shown as filled circles in Fig. |l^(b), 
exhibits good agreement with the asymptotic law. As for the values discussed for the tagged- 
particle density correlators, the drift of v^^ still results in a temperature dependence of /l^^e 
for T > Tc, but this drift is now reduced to a 10% effect. The asymptotic value of f^^^^, dif- 
fers only about 5% from the one for the first correlator, f^, which is shown as a dash-dotted 
line in Fig. |T^(b). It is remarkable, that even for the unnormalized quantity the position of 
Tc can be estimated better than it could be done for the plateau values of the tagged-particle 
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density correlators. This can be done by noticing that the slope of a linear interpolation of 
the data changes when going over from T < to T > T^. 

From the MCT equations (|1]) with the hopping kernel set to zero, one derives the ex- 
pression for the dynamic susceptibility, Eq. (|^), in terms of the memory kernel m^(z). 



Xa{z) = -^WaI 



(21) 



Let us define a dynamical susceptibility Xm,Ac{A corresponding to the kernel rn%Sf) in 
analogy to Eq. (|^): 



(22) 



with m^g = "^de(^ — 0)- Then one can write for the dielectric function e{z^ = £oo+4vrxde(^) 

-47rxde 



Biz 



(23a) 



(z/fiSJ^-l-m^,,o + 

The inverse of the dielectric function, l/e{z), is occasionally considered as the dielectric 
modulus [Q. The exact Mori-Zwanzig representation, Eq. (0), suggests to rather con- 
sider [e{z) — ^oo]"^! i- e. Xdc(-2)- This function consists of a quadratic polynomial in the 
frequency, (^/fi^^)^ — 1 — m^gQ, a white-noise background, izu^^, and a non-trivial part 
XM,Ae{z). The latter has all the standard properties of a susceptibility, in particular it obeys 
Kramers-Kronig relations. There is the trivial relation between the spectrum xX/dol"^) ^^"^ 
the dielectric function 



Im 



1 



4vrxdc 



del + 



(23b) 



Here, £00 was replaced by the constant e, discussed above in connection with the fit of s'{uj). 

The full lines in Fig. ^ exhibit the right-hand side of Eq. ( |23 b| ) , evaluated with the 
model parameters used for the interpretation of the dielectric-loss spectra in Fig. |1|. The 
symbols exhibit the left-hand side of Eq. ( p3b| ) calculated with the data from Ref. and 
e determined in connection with the fits of s'{uj) in Fig. ^ Figure |TB| shows that the fit 
is of equal quality as the ones shown for the direct analysis of the dielectric loss spectra. 
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However, to produce the result in Fig. |T^, one has to be careful to subtract the right value 
of i. The error bars shown in the figure for T = 253 K indicate the influence of subtracting 
e ± 1 instead of e to estimate the uncertainty introduced by this procedure. One notices that 
the shape of the curves for the high-frequency part is influenced. Thus an analysis based on 
[e — e{uj)]~'^ is only practicable, if one can avoid these problems of the inversion procedure. 
Up to trivial terms, the left-hand side of Eq. (|23b| ) is identical with the spectrum um'^^^u) 
discussed in Fig. This quantity can be explained well by the /3-relaxation-scaling laws. 
Thus for PC, the corrections to the asymptotic laws are smaller for ( [£(^1;) — i] 1 than for 



s"{uj)- In particular it is shown by Fig. |T5|that the square-root singularity can be identified 
from a discussion of XM,de{^^)- 



VII. CONCLUSION 

It was exemplified for propylene carbonate as a typical glass-forming van-der-Waals sys- 
tem, that the susceptibility spectra measured by three different experimental techniques can 
be described well by a schematic MCT model. Several decades of intensity change in the 
GHz frequency window as seen in light-scattering and dielectric-loss experiments with a sen- 
sitive temperature dependence typical for glass-forming liquids are fitted. Also, the results 
from incoherent neutron scattering, probing the dynamics for different wave vectors could be 
included in this simultaneous fit. Real-part data from the dielectric experiment have been 
successfully analyzed as well to further corroborate the consistency of the schematic-model 
fit. For temperatures ranging from the critical value Tc ~ 180 K to well above the melting 
point, the range of applicability of the model includes both a- and /^-relaxation windows, as 
well as the crossover to the microscopic spectrum. Below Tc and down to the glass-transition 
temperature Tg, a rather simple approach to account for hopping phenomena improves the 
fit for the /3-minimum regime, but fails to describe the a-peak below Tc. 

The schematic model used in the fit captures the general features of the glass-transition 
scenario predicted by the full microscopic theory. Still, it allows to go further than an anal- 
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ysis based on the asymptotic predictions of MCT only. In particular, we have used the 
schematic model to investigate which features of the measured spectra can be described by 
asymptotic laws, and where preasymptotic corrections set in. We find that the asymptotic 
formulas qualitatively give an adequate description of the data. Thereby the preceding stud- 
ies P],|3|-|7|,pi| are corroborated. But we demonstrated also, that preasymptotic-correction 
effects cause important quantitative differences between the data and the scaling-law results. 
One aspect of this is the T-drift of the critical amplitude hf^ noted in an earlier analysis 
of PC light-scattering data The drift of the coupling constant ffg is not sufficient to 
explain this, as was demonstrated in Figs. ^ and |T3|. Also, the crossover to the microscopic 
excitations influences the height of the spectra at the /3 minimum. For the measurements 
analyzed, scaling works best with the neutron-scattering data, due to the relatively low 
plateau values f^'^. An asymptotic analysis of the dielectric modulus could even work better 
in this respect. But due to uncertainties in the inversion of the dielectric function, such 
analysis is not practicable unless the modulus itself is measured directly. 



ACKNOWLEDGMENTS 

We thank H. Z. Cummins, M. Fuchs, P. Lunkenheimer, M. R. Mayr, U. Schneider, A. P. 
Singh, and J. Wuttke for many helpful discussions and the authors of Ref. [Q, Refs. 
and Ref. [|l^] for providing us with their files for the various propylene-carbonate data. This 
work was supported by Verbundprojekt BMBF 03-G05TUM. 



30 



REFERENCES 

[1] W. M. Du, G. Li, H. Z. Cummins, M. Fuchs, J. Toulouse, and L. A. Knauss, Phys. Rev. 
E 49, 2192 (1994). 

[2] W. Gotze and L. Sjogren, Rep. Prog. Phys. 55, 241 (1992). 

[3] J. Ma, D. V. Bout, and M. Berg, Phys. Rev. E 54, 2786 (1996). 

[4] P. Lunkenheimer, A. Pimenov, M. Dressel, B. Gorshunov, U. Schneider, B. Schiener, 
and A. Loidl, Am. Chem. Soc. Symp. Ser. 676, 168 (1997). 

[5] P. Lunkenheimer, A. Pimenov, M. Dressel, B. Schiender, U. Schneider, and A. Loidl, 
Prog. Theor. Phys. Suppl. 126, 123 (1997). 

[6] U. Schneider, P. Lunkenheimer, R. Brand, and A. Loidl, Phys. Rev. E 59, 6924 (1999). 

[7] L. Borjesson, M. Elmroth, and L. M. Torell, Chem. Phys. 149, 209 (1990). 

[8] M. Elmroth, L. Borjesson, and L. M. Torell, Phys. Rev. Lett. 68, 79 (1992). 

[9] W. Gotze, J. Phys.: Cond. Matt. 11, Al (1999). 

[10] J. Wuttke, M. Ohl, M. Goldammer, S. Roth, U. Schneider, P. Lunkenheimer, R. Kahn, 
B. Ruffle, and R. Lechner, Phys. Rev. E, 1999, submitted. 

[11] N. V. Surovtsev, J. A. H. Wiedersich, N. V. Novikov, E. Rossler, and A. P. Sokolov, 
Phys. Rev. B 58, 14888 (1998). 

[12] H. C. Barshiha, G. Li, G. Q. Shen, and H. Z. Cummins, Phys. Rev. E 59, 5625 (1999). 

[13] T. Pranosch, M. Fuchs, W. Gotze, M. R. Mayr, and A. P. Singh, Phys. Rev. E 55, 7153 
(1997). 

[14] M. Fuchs, W. Gotze, and M. Mayr, Phys. Rev. E 58, 3384 (1998). 

[15] U. Bengtzelius, W. Gotze, and A. Sjolander, J. Phys. C 17, 5915 (1984). 



31 



[16] W. Gotze and L. Sjogren, Z. Phys. B 65, 415 (1987). 

[17] F. Sciortino and P. Tartaglia, J. Phys.: Cond. Matt. 11, A261 (1999). 

[18] C. Alba-Simionesco and M. Krauzman, J. Chem. Phys. 102, 6574 (1995). 

[19] C. Alba-Simionesco, V. Krakoviack, M. Krauzman, P. Mighardo, and F. Romain, J. 
Ram. Spect. 27, 715 (1996). 

[20] V. Krakoviack, C. Alba-Simionesco, and M. Krauzman, J. Chem. Phys. 107, 3417 
(1997). 

[21] W. Gotze, Z. Phys. B 56, 139 (1984). 
[22] L. Sjogren, Phys. Rev. A 33, 1254 (1986). 

[23] T. Franosch, M. Fuchs, W. Gotze, M. R. Mayr, and A. P. Singh, Phys. Rev. E 56, 5659 
(1997). 

[24] T. Franosch, W. Gotze, M. R. Mayr, and A. P. Singh, Phys. Rev. E 55, 3183 (1997). 

[25] A. P. Singh, G. Li, W. Gotze, M. Fuchs, T. Franosch, and H. Z. Cummins, J. Non-Cryst. 
Solids 235-237, 66 (1998). 

[26] B. Ruffle, C. Ecolivet, and B. Toudic, Europhys. Lett. 45, 591 (1999). 

[27] W. Gotze, in Liquids, Freezing and Glass Transition, edited by J. P. Hansen, D. 
Levesque, and J. Zinn- Justin (North-Holland, Amsterdam, 1991), pp. 287-503. 

[28] M. Fuchs, W. Gotze, S. Hildebrand, and A. Latz, J. Phys.: Cond. Matt. 4, 7709 (1992). 

[29] T. Franosch, W. Gotze, M. R. Mayr, and A. P. Singh, J. Non-Cryst. Sohds 235-237, 
71 (1998). 

[30] P. B. Macedo, C. T. Moynihan, and R. Bose, Phys. Chem. Glasses 13, 171 (1972). 



32 



FIGURES 

FIG. 1. Susceptibility spectra for propylene carbonate (PC, ~ 218 K, Tg = 160 K) as 
measured by depolarized-light-scattering (upper panel, data from Ref. Q), and dielectric- loss 
spectroscopy (lower panel, data from Ref. Q), normalized with a temperature- independent static 
susceptibility. Temperatures are in steps of 10 K, unless indicated otherwise; in the dielectric mea- 
surement, T = 243, 263, 273, 283 K are missing, for the light-scattering experiment, the highest 
temperature is 250 K (see text for details). The full lines are fits by solutions of the two-component 
schematic MCT model defined in Sec. || with parameters as described in the text. The dashed lines 
indicate a white-noise spectrum, n i^) ^ '^i dash-dotted line in the upper panel exhibits the 
asymptote of the critical spectrum ^"{u)) oc according to Eq. ([l5|), with a = 0.30 corresponding 
to the value of A = 0.75. The dotted line shows the solution of the model for T = Tc and hopping 
terms neglected. 

FIG. 2. Susceptibility spectra for PC as measured by incoherent neutron scattering for three 
wave vectors q, from Ref. 0. Temperatures are T = 210, 220, 230, 240, 251, 260, 285 K, where 
alternating open and filled symbols have been used to help distinguishing different data sets. 
For lower temperatures, data points below 10 GHz are seriously affected by the spectrometer's 
resolution function and therefore not shown. Full lines are fits as in Fig. ||. 

FIG. 3. Vertices vi, V2 for the first mode-coupling functional, Eq. (^), used for the fits shown 
in Figs. 1^ and ^. In the above diagram, the thick line represents the curve of glass-transition 
singularities, Eq. (^), while the thin line serves as a guide to the eye indicating the chosen path. 
Each dot corresponds to one temperature. The lower diagram shows the separation parameter a, 
Eq. (P), as a function of T; the critical temperature ~ 180 K is determined from the zero of the 
shown regression line. 
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FIG. 4. (a) Coupling coefficients v% for the second mode-coupling functional, Eq. (^), used for 
the fits shown in Figs. Q and ^ as functions of temperature. The squares refer to the light-scattering 
data and the circles to the dielectric-loss spectra; the lines through the symbols are guides to the 
eye. The lines without symbols connect the v^^{q) used for the neutron-scattering data for 10 wave 
vectors q = 0.5, 0.6, 1.4 A ^ (from top to bottom). The vertical dashed line indicates the 
critical temperature Tc. (b) Coefficients v^^{q,T) for the neutron-scattering data as functions of q 
for various fixed T. From top to bottom, the temperatures increase from T = 210 K to 285 K (as 
given in Fig. |2|). The dashed lines indicate A/ q + B / q^ laws to visualize the difference to a 1/g-law 
behavior, which is shown as a dot-dashed line; see text for details. 

FIG. 5. Oscillator frequencies and damping constants v\ in THz specifying the transient 
motion for the second correlator as functions of temperature as used for the fits shown in Figs. |l| 
and|2|. The squares refer to the light-scattering spectra, circles to the dielectric loss spectra; lines 
through the symbols are guides to the eye. For the neutron-scattering spectra, z^*g was taken 
^-independent (diamonds in the lower panel), and the fi^g exhibit the q ■ v^-behavior shown in 
the upper panel (lines without symbols, q range as in Fig. |^. The vertical dashed lines indicate 
the critical temperature Tc. 

FIG. 6. In the upper part, the hopping coefficient A entering Eq. ( |3a[ ) for the ffi'st correlator 
used for the fits shown in Fig. |^ to the dielectric (circles) and light-scattering (squares) data is 
exhibited. The values follow an Arrhenius-type temperature dependence, indicated by a straight 
line. For the fit to the neutron-scattering data, the same values could be used, but show no influence 
on the fit curves in Fig. |2| (see text for details). For the fit to the dielectric data, an additional 
hopping coefficient A^^ for the second correlator had to be used, shown by the circles in the lower 
part of this figure; here, the values (J7jp/r2)^A^p follow an Arrhenius law indicated by a straight 
line. The vertical dashed line indicates 1000/Tc with Tc = 180 K. 
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FIG. 7. The solid lines reproduce the susceptibility spectra x"a{'^) from Fig. [l| used for the 
fit to the light-scattering (upper panel) and the dielectric- loss spectra (lower panel), respectively. 
Dashed lines are solutions using the same model parameters, but with hopping effects ignored: 
A = = 0. Notice that the dashed lines for temperatures T below Tc exhibit a "knee" which is 
located between 10 and 100 GHz and moves to higher frequencies with decreasing T. 

FIG. 8. Measured values for the real part e'{uj) of the dielectric function from Ref. |^] for 
T = 173 K through T = 233 K in steps of 10 K, and for T = 253 K (from left to right). Exper- 
imental data have been shifted by e to account for an unknown background; see text for details. 
Temperatures T = 163 K and below, and T = 293 K are not shown in order to avoid overcrowding 
of the figure. For the same reason, the data points for frequencies above 1 GHz have been partially 
removed for all but the highest temperature and are only shown in the inset. The full lines are the 
real parts of the calculated susceptibilities for the same model parameters as used for the curves 
in Fig. ||. 

FIG. 9. The full lines show 4)x{t) = i4>xit) — fx)/hx for times where cpxit) > 0. The dash-dotted 
straight lines exhibit (t/to)~") with Iq = 0.035 ps and a = 0.30. The dashed lines show the 
/3-correlators G{t) for the indicated temperatures (see text). The uppermost set of curves refers 
to the first correlator of the schematic model. The other two sets refer to the light-scattering and 
the dielectric response studied in Fig. ||. The sets are shifted vertically by one (light scattering), 
respectively two (dielectric response) decades for clarity. 
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FIG. 10. Effective nonergodicity parameters determined from the calculated (j)x-veTsus-logt 
curves for the spectra shown in Fig. ^ f{T) for the first correlator (j){t) (open circles), and f%{T) 
for the correlators 0fg(i) and 0jg(t) (filled squares and circles, respectively). The solid lines exhibit 
the leading-order-asymptotic laws for fx from Eq. (^) for T < Tc and for T > T^. The vertical 
dashed line indicates the crossover temperature Tc = 180 K. The crosses in the lowest panel are 
estimations of f^^ obtained in Ref. by fitting the measured s'{uj) data with a Cole-Davidson 
function. The dashed lines through the solid symbols connect the values calculated from Eq. ( |18bD 
(see text for details). 

FIG. 11. The upper part of the figure (scale on right axis) shows the rescaled susceptibility 
spectra for four temperatures for the first correlator of the schematic model, uj(p" {uj)/h^/\a\, as 
a function of cu = uji^j. Hopping terms have been set to zero. The lower part (scale on the 
left axis) shows the tiJ(/)^"(w)//iA\/|o^-versus-cj curves for T = 210 K for the light-scattering (line 
with squares) and neutron-scattering correlators at q = 1.0 A ^ and q = l.lA ^ (lines), and for 
T = 213 K for the dielectric correlators (line with circles). The dashed lines indicate the master 
spectrum for A = 0.75 with the minimum position (a)min, Xmin) marked by a diamond. 

FIG. 12. Rectification diagram for the minimum positions Lo^iin of the dielectric-loss (circles), 
light-scattering (squares), and neutron-scattering (triangles) fit. The neutron-scattering data have 
been averaged over the investigated q range, and error bars indicate the smallest and highest values. 
The vertical dashed line marks the critical point Tc- The solid line is the asymptotic-formula result 
for the minimum positions calculated from the /3-relaxation-scaling law for the fit parameters of the 
model used in Figs. ^ and|^. The dashed lines show the linear interpolation through the data points, 
restricted to the regime 190 K <T< 250 K, for the dielectric-loss and for the light-scattering data, 
respectively. 
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FIG. 13. Scaling-law analysis of the MOT solutions shown in Fig. |l| for the fit of the 
light-scattering spectra. The lower set of curves shows the rescaled spectra X\s(.^)/^isVW\ fo^' 
T = 220, 210, 200, 190 K (from top to bottom). The upper set of curves shows the same spectra, 
but now rescaled with his replaced by the effective amplitude hf^ from Eq. (|18b| ). Notice that the 
order of the rescaled spectra is now inverted. The asymptotic result x is plotted as a dashed line 
with the minimum position marked by a diamond. The result for T = 180 K with hopping effects 
ignored is added to demonstrate the approach towards the scaling limit. 

FIG. 14. Scaling-law analysis of the spectra for the memory kernel from the MOT solutions 
shown in Fig. |^ to interpret the dielectric function. The upper set of solid curves shows the 
unnormalized rescaled spectra, —u;m'^^{uj) / ^/\a\, for temperatures T = 193, 203, 213 K. In the 
lower part of the figure, the scaling is shown for the normalized rescaled spectra —ujm'^^{uj) / y^Ja\ 
(see text for details). The dashed lines are the master spectra, with the minimum positions marked 
by diamonds. 

FIG. 15. (a) Unnormalized nonergodicity parameters fl^^^ for the memory kernel from the 
fit to the dielectric spectra (open circles). The asymptotic prediction is plotted as a solid line, the 
vertical dashed line indicates Tc. (b) Same as in (a), but normalized values /|,/(je (fiH^d circles). 
The dash-dotted line indicates the critical plateau value for the first correlator of the schematic 
model, f^. 

FIG. 16. The lines show the dielectric-modulus spectra defined by the right-hand side of 
Eq. (23b). The circles exhibit Im [e — e(Lj)]^^ 



calculated from the dielectric-function data of 



Ref. 1^. The curves and the data sets have been shifted vertically to avoid overlapping; tempera- 
tures and shift factors are (183 K, 1), (193 K, 0.316), (203 K, 0.1), (213 K, 0.0316), (223 K, 0.01), 
(233 K, 0.00316), (253 K, 0.001). For the highest temperature, an estimation of the uncertainty 
introduced by inverting the experimental data is given as error bars; see text for details. 
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Gotze and Voigtmann Figure 1 
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